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We compute cross sections for incoherent diffractive J/'i' production in lepton- nucleus deep in- 
elastic scattering. The cross section is proportional to A in the dilute limit and to A^''^ in the black 
disc limit, with a large nuclear suppression due to saturation effects. The f-dependence of the cross 
section, if it can be measured accurately enough, is sensitive to the impact parameter profile of 
the gluons in the nucleus and their fluctuations, a quantity that determines the initial conditions 
of a relativistic heavy ion collision. The nuclear suppression in incoherent diffraction shows how 
the transverse spatial distribution of the gluons in the nucleus gradually becomes smoother at high 
energy. Since the values of the momentum transfer |t| involved are relatively large, this process 
should be easier to measure in future nuclear DIS experiments than coherent diffraction. 
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I. INTRODUCTION 

Strongly interacting systems in the high energy (or 
small x) limit are very nonlinear systems in spite of the 
smallness of the coupling constant . This is due to the 
large phase space available for semihard gluon radiation 
that increases the occupation numbers of gluonic modes 
in the hadron or nucleus wavefunction. Thus high en- 
ergy scattering has to be understood in terms of gluon 
recombination and saturation that enforce the unitarity 
requirements of the ^-matrix. This happens naturally in 
the Color Glass Condensate (CGC) effective theory of the 
high energy wavefunction. In the context of deep inelas- 
tic scattering (DIS) the CGC leads to the dipole picture 
that naturally gives a consistent description of both in- 
clusive and diffractive scattering. The nonlinearities in 
high energy scattering are enhanced when the target is 
changed from a proton to a heavy nucleus. Thus there is 
a great opportunity to understand them by studying nu- 
clear DIS in new collider experiments, such as the EIC [T] 
or the LHeC [5] . The particular process we discuss in this 
paper is diffractive DIS on nuclei. 

In the Good- Walker |3] picture of diffraction one needs 
to identify the states that diagonalize the imaginary part 
of the T-matrix. In the case of nuclear DIS at high en- 
ergy these states are the ones with the virtual photon 
fluctuating into a dipole of a fixed size r and with the 
nucleons in the nucleus at fixed transverse positions bi. 
In coherent diffraction the nucleus is required to stay in- 
tact, which corresponds to performing the average over 
the nuclear wavefunction at the level of the scattering 
amplitude. Averaging the cross section, instead of the 
amplitude, over the nucleon positions allows for the nu- 
cleus to break up, giving the sum of incoherent and co- 
herent cross sections, i.e. the quasielastic cross section. 
For a more formal discussion of this we point the reader 
e.g. to Ref. [5. The i-dependence of the incoherent cross 
section therefore directly probes the fiuctuations and cor- 
relations in the nuclear wavefunction, which have turned 
out to be a crucial ingredient in understanding the initial 



conditions in heavy ion collisions [S]. 

The average gluon density probed in the coherent pro- 
cess is very smooth, meaning that the cross section is 
dominated by small values of momentum transfer to the 
nucleus, t ^ —1/R\. Measuring such a small momen- 
tum transfer accurately is very challenging. At momen- 
tum scales corresponding to the nucleon size t ^ — 
the diffractive cross section is almost purely incoherent. 
The larger momentum transfer should also be easier to 
reconstruct experimentally even without measuring the 
transverse momentum of the nuclear remnants, by ac- 
curately reconstructing the outgoing electron and J/^ 
momenta and using momentum conservation. By taking 
these processes into account in the detector design one 
should be capable of measuring diffractive events at a 
higher accuracy than was done at HERA. In the dilute 
limit (for small dipoles) there is no multiple scattering, 
and the incoherent cross section is given by A times the 
corresponding one for protons. The deviation of the t- 
slope from the proton measures the transverse size of the 
fluctuating areas in the nucleus. 

In the black disc limit the nucleus is smooth not only on 
average, but event-by-event, leading to a strong suppres- 
sion of the incoherent cross section. Incoherent diffrac- 
tion gets contributions from the edge of the nucleus, mak- 
ing the cross section asymptotically behave as ^ A^^^ in 
contrast to ~ A in the dilute limit. The suppression in 
the normalization relative to the proton is a measure of 
the approach to the unitarity limit in the dipole cross 
section. It is a clear signal of how individual nucleons 
have lost their identity in the sense that they cannot be 
resolved by the virtual photon. It is precisely this sup- 
pression that we are proposing to use to quantitatively 
access saturation effects in the nuclear wavefunction. The 
purpose of this paper is to provide a realistic estimate of 
the nuclear suppression in diffractive cross sections in a 
regime that could be measured in future nuclear DIS ex- 
periments. 

Nuclear DIS data from fixed target experiments, in 
particular E665 [6] and NMC [7] have already been much 
discussed in the literature as demonstrations of color 
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transparency (see e.g. Refs. [8til2)V The form of nu- 
clear modification to the incoherent diffraction in terms 
of the dipole cross section that we have rederived is not 
new (see e.g. [HI US])- So far, however, less attention has 
been paid to inelastic diffraction in future DIS experi- 
ments. The production cross sections have not been cal- 
culated using the same CGC inspired cross sections that 
have been used successfully to confront HERA data, as 
we intend to do here. In this work we concentrate on the 
J/vl/ because its small size means that the interaction of 
the dipole with the target is calculable in weak coupling 
even at small Q^. 

The importance of diffraction in understanding gluon 
saturation has been discussed and our basic setup moti- 
vated in Ref. [IT. Nuclear modifications to the diflFrac- 
tive structure functions, integrated over the momentum 
transfer t, were computed in Ref. [15 . Vector meson 
production at future DIS experiments was recently dis- 
cussed from a more experimental point of view in Ref. 4J , 
and coherent production cross sections (integrated over 
t) calculated in Ref. [16]. An interesting discussion on 
coherent and incoherent diffraction and gluon saturation 
in the nucleus can be found in Ref. |17j . In this study 
we want to take a step beyond the discussion of inclusive 
diffraction in Refs. [Ill [15] to understand the t depen- 
dence in more detail. 



II. DIPOLE CROSS SECTIONS 

There are many dipole cross section parametrizations 
available in the literature, and we have taken for this 
study two representative samples. One is the IIM [TS] 
dipole cross section, which is a parametrization includ- 
ing the most important features of BK [T^] evolution. 
The detailed expression for the dipole cross section can 
be found in Ref. [TH]; we use here the values of the pa- 
rameters from the newer fit to HERA data including 
charm [5D] that was also used to compute diffractive 
structure functions in Ref. 21 . Wc also want to com- 
pare our results to a parametrization with an eikonalized 
DGLAP-evolved gluon distribution. For this purpose we 
will use an approximation of the IPsat dipole cross sec- 
tion [ll[23j. 

To extend the dipole cross section from protons to nu- 
clei we will take the independent scattering approxima- 
tion that is usually used in Glauber theory and write the 
S'-matrix as 
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(1) 



Here we conventionally parametrize the energy depen- 
dence of the scattering amplitude with x, the Bjorken 
lable of the DIS event^. The variables hn in Eq. ^ 
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FIG. 1: The r-dependence of the different proton dipole 
cross sections used, at a; = 0.0001 and & = 0. As dis- 
cussed in Sec. IV the "IPnonsat" -curve is Eq. (|5| linearized 



^F{x,r). 



are the nucleon coordinates that we will discuss in 
Sec.|III| This independent scattering assumption is natu- 
ral in IPsat-like parametrizations or the MV [23] model, 
where, denoting r = Ir^l, S{rT) ~ "^^^^ with a sat- 
uration scale Qg proportional to the nuclear thickness 
TA{b). High energy evolution, however, introduces an 
anomalous dimension that leads, in the nuclear case, to 
what could be called leading twist shadowing. With an 
anomalous dimension S ~ e"'''^^''-' ^ with 7 7^ 1, a pro- 
portionality ~ TA{b) is not equivalent to Eq. ([l]). 
A solution to this problem (see also the more detailed 
discussion in |15j ) would require a realistic impact pa- 
rameter dependent solution to the BK equation which, 
we feel fair to say, is not yet available. We point the 
reader e.g. to Ref. [3S] for a discussion of the difficulties. 
These are related to the long distance Coulomb tails that, 
physically, are regulated at the confinement length scale 
that is not enforced in a first principles weak coupling 
calculation. The effect of BK evolution is important for 
the CGC description of the forward suppression of par- 
ticle production in dAu-collisions at RHIC (for a review 
see [26j). In our case the difficulty is greater since we 
are interested not only in the relatively smooth average 
gluon density, but its variations at smaller length scales 
of the order of the proton radius. We thus leave the mod- 
ifications of Eq. ([ij due to the effects of evolution to a 
future study. 

The IIM parametrization assumes, either explicitly or 



Note that strictly speaking the relation between x and the en- 



ergy of the dipole-target scattering depends on , not only r. 
Using X here is justified in a high energy approximation where 
the energy of the dipole in the target rest frame is approximately 
the same as that of the virtual photon. 



3 



implicitly, a factorizable by dependence 



7*p^,J/*p, W=90 GeV 



d^hr 



(bT,rT,x) = 2(l-5p(rT,bT,a;)) (2) 
= 2rp(bT)AA(r,a;), 



We take, following Ref. [2T], a Gaussian profile TJhx) = 
exp {-b^/2Bp) with Bp = 5.59 GeV-^ (see Sec.^for a 
discussion of this largish numerical value). 

In the IPsat model the impact parameter dependence 
is included in the saturation scale as 



da 



d^h 



^{hr, TT, x) = 2 [1 - exp {^r'F{x, r)Tp(bT))] 



(3) 

Here TpCbx) = exp (— &^/2i3p) is the impact parameter 
profile function in the proton with Bp = 4.0 GeV^ and 
F is proportional to the DGLAP evolved gluon distribu- 
tion [57] 
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FIG. 2: Comparison of the used dipole cross sections to HERA 
data [28] on diffractive vector meson production. 



III. COMPUTING DIFFRACTIVE CROSS 
SECTIONS 
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(4) 

with C chosen as 4 and /i§ = 1.17 GeV resulting from 
the fit [22]. The proton dipole cross sections used are 
plotted in Fig.[l]for x = 0.0001. 

We would generally prefer the unfactorized b- 
dependence of Eq. ([3| to the factorized one in Eq. ^ 
because it allows for the correct unitarity limit of the 
scattering amplitude at all impact parameters (see the 
discussion in Ref. [15]). However, there seems to be no 
clear difference between the two in terms of the quality of 
the description of HERA data, and for the sake of com- 
putational simplicity we will in this work limit ourselves 
to the factorized dependence and approximate the IPsat 
dipole cross section by 



d2bT 



{hT,rT,x) 



2Tp(bT) [l-exp{-r''F(x,r 



))] 



(5) 

using the same F{x,r) defined in Eq. Q. This approx- 
imation brings the IPsat parametrization to the form 
Eq. |2| with7V(r, x) = [l - exp (-r2i^(a;, r))] ; in fact 
this IS the form used already in Ref. |27| : we however 
use the gluon distribution from the IPsat fit [23] for con- 
venience. Improving this description goes hand in hand 
with giving up the approximation of independent scat- 
terings off the nucleons, Eq. ([T]), and is left for future 
work. As we shall see in the following, these approxima- 
tions enable us to write the cross section for incoherent 
diffraction in a form which is much simpler to evaluate 
numerically than one with a general 6-dependence. 



The cross section for quasielastic vector meson produc- 
tion in nuclear DIS is 



dCT^*^- 



■VA 
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(6) 



with t = —At ■ The dipole cross section is evaluated 
at the energy scale corresponding to the rapidity gap be- 
tween the vector meson and the target xp. To translate 
this into the photon-target center of mass energy W that 
is often used to present experimental results note that 
xp = (Mj/^ + Q^)/{W^ + Q2). The factor 1 + (3^ ac- 
counts for the real part of the scattering amplitude and 
the factor corrects for the skewedness effect, i.e. that 
the gluons in the target are probed at slightly different 
X ^29j . For these corrections we follow the prescription of 
Ref. [30] J taking them as 



= tan 



ttA 



Ro 



A + 4 



with 



A = 



>J/'i'P 



d In 1 /xp 



(7) 
(8) 

(9) 



These corrections depend, in general, on t, which we 
take into account in our calculation. For the full IP- 
sat model A changes by about 5% between t = and 
— t = 0.5 GeV^. For the factorized impact parameter de- 
pendence in Eqs. ^ and ([5]) A is independent of t. We 
calculate the correction terms from the energy depen- 
dence of the nucleon scattering amplitudes and use the 
same values for the nucleus at the same Q'^^xp. Since 
the difference in A extracted from the nucleus and the 
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nucleon cross sections is small (compared to the value of 
A) and Rg and /? are in themselves corrections to the cross 
section, this approximation is justified. In addition this 
approximation has the advantage that these corrections 
cancel on the nucleus/nucleon cross section ratio. The 
real part and skewedness corrections, especially Rg are, 
however, a significant factor in the absolute normaliza- 
tion of the cross section and are necessary for the agree- 
ment with HERA data. 

The imaginary part of the scattering amplitude is the 
Fourier-transform of the dipole cross section from hx to 
A.T contracted with the overlap between the vector me- 
son and virtual photon wave functions: 



A = 197, Q =0 GeV^ a;p=0.001 



X ['i>m{r,Q',z)e~^^- '^-^{hT,rT,xr), (10) 

where we have followed the normalization convention 
of |23j . For the virtual photon-vector meson wavefunc- 
tion overlap we use the "boosted Gaussian" parametriza- 
tion from Ref. [25]. We have also tested the "gaus- 
LC" wavefunction also used in Ref. [53]. Although the 
"boosted Gaussian" seems preferred by HERA data, also 
the "gaus-LC" parametrization is compatible with the 
data within the experimental errors. The cross sections 
for the proton differ by factors of the order of 10%. The 
interaction of the gluon target with the dipole can in 
general depend also on which introduces terms that 
couple Ft, At and z in Eq. (10). For the J/^E* and the 



range in t considered in this paper Ay is sufficiently small 
compared to the relevant values of 1/r that we can ne- 
glect this coupling, which simplifies the structure con- 
siderably. Lighter vector mesons would require a more 
general treatment. 

The average over the positions of the nucleon in the 
nucleus is denoted here by 



(0({bTj))r 



n[d'bT,;T^(bT.)]0({bTj). (11) 



Here Ta is the Woods-Saxon distribution with nuclear ra- 
dius Ra = (1.12^1/3 _ 0.86A-1/3) fm and surface thick- 
ness d = 0.54 fm. This expectation value is equivalent to 
the average over nucleon configurations in a Monte Carlo 
Glauber calculation. We are assuming that the positions 
b^i are independent, i.e. neglecting nuclear correlations 
that would be a subject of interest in their own right (see 
e.g. [3T]). The coherent cross section is obtained by av- 
eraging the amplitude before squaring it, | (^)i,^ and 



the incoherent one is the variance {\A\'^ 



(A), 



that 



measures the fluctuations of the gluon density inside the 
nucleus. Because {A}-^ is a very smooth function of by, 
its Fourier transform vanishes rapidly for A > 1/Ra- 
Therefore at large A the quasielastic cross section ^ is 
almost purely incoherent. 

The cross section for quasielastic vector meson produc- 
tion is now expressed in terms of the dipole scattering 




Quasielastic IPnonsat 
Colnerent IPnonsat 
IPsat 

Coherent IPsat 
MM 

Colnerent MM 



10^ 



0.05 



0.10 



0.15 



0.20 



0.25 



|i| [GeV' 



FIG. 3: The quasielastic and coherent diffractive J/^ cross 
sections in gold nuclei at = and xp = 0.001. Shown are 
the IPsat and IIM parametrizations. We also show the result 
for the linearized "IPnonsat" version (used e.g. in Ref. [3]) 
where the incoherent cross section is explicitly A times that 
of the proton. Our approximation ( |13[ ) is not valid for small 
\t\; the corresponding part of the distribution has been left 
out. 



amplitude as 



dt 167r J All Att 

X [^*y^ (r,z,Q) [^*y^ (/,z',Q) 

lA?!^ {xv,r, r', At; 



d^rT d^rT' 



(12) 



We now average the square of the dipole scattering am- 
plitude over the nucleon coordinates, using the assump- 
tions of Eqs. ([!]) and ^ and taking the large A limit. 
We are additionally assuming that T4 is a smooth func- 
tion on the discance scale defined by Bp. Averaging the 
square of the amplitude gives the total quasielastic con- 
tribution, but we only keep the terms that contribute at 
large \t\ ^ 1/-R^, which leaves us with the expression 



A 

\Agii\' (xp, r, r'. At) = 167rBp j d^hr ^ 



n \n 

^ g-BpATV"e-2'rSpATA(6)[AA(r)+7V(r')] 

nBpM{r)M{r')TA{b) 



l-2^BpTA{b)W(r)+N{r')] 



(13) 



Note that Eqs. ([!]) and ^ have enabled us to write 
the leading contributions as proportional to the (Gaus- 
sian) proton impact parameter profile, which can then 
be Fourier-transformed analytically. Giving up either 
of these approximations would force us to numeri- 
cally Fourier-transform the "lumpy" 6-dependence cor- 
responding to a fixed configuration of the nucleon posi- 
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tions. This would make the numerical calculation much 
more demanding and is left for future work. 

The terms with n > 2 correspond to scattering off a 
system of several overlapping nucleons simultaneously, 
leading to slower suppresion with \t\. In practice we have 
verified numerically that they do not contribute to our 
results at the values of t we are interested in (the n — 2 
contribution is typically < 2% of the n — 1-one, only 
reaching 5% at —t > 0.5 GeV^ ) and will neglect them 
in the following. This leaves us with the expression 

\Aqgfixp,r,r',AT)^WTTBpA J d^hr 

^ ^-BpA.T\-27rBpATA(b)[\rir)+Af{r')] 



A =197, \t\ =0.5 GeV^ 



7TBpJ^{r)N{r')TA{b) 
l-27rBpTA{b)[J\fir)+J\f{r')] 



(14) 



Equation ( 14 1 has a very clear interpretation. The 
squared amplitude is proportional to A times the squared 
amplitude for scattering off a proton, corresponding to 
the dipole scattering independently off the nucleons in 
a nucleus. This sum of independent scatterings is then 
multiplied by a nuclear attenuation factor 

^~^27rB^ATA{b)[Af{r)+Af(r')] 

l-2iTBpTA{b)[Af{r)+Af(r')] ^ 

^-2wiA-l)B^TA{b)[N-{r)+N-(r')] ^ ^--^g-j 

which accounts for the requirement that the dipole must 
not scatter inelastically off the other A — 1 nucleons in 
the target (otherwise the interaction would not be diffrac- 
tive). Note that the factor 4TTBpAf{r,xp) = cr^jp(r, xp) is 
the proton-dipole cross section for a dipole of size r. Thus 
this attenuation corresponds to the probability of a dipole 
with a cross section which is the average of dipoles with 
r and r' to pass though the nucleus. A similar expression 
can be found e.g. in Ref. [TTj . 

For comparison, the coherent cross section in our ap- 
proximation is given by 



dt 



167r 



{A{xp,Q\At))J\ (16) 



where in the large A and smooth nucleus limit the am- 
plitude is 
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FIG. 4: The "nuclear transparency" ratio of cross sections 
vs. for IPsat, IIM parametrizations at xv = 10~^ (the 
upper three curves, blue) and 10"'' (the lower 3 curves, black). 
For comparison we also include we also include the result if 
unitarization effects are included at the nucleus but not at 
the nucleon level in the IPsat-parametrization. (See text for 
discussion). 
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using the "Gaus-LC" vector meson wavefunctions. The 
labeling is the same as in Fig. [4] 



{A{xpM\^t))^^ I ^dVd^bTC-b-^- 

X [^*y^]ir, g2, z) 2 [1 - exp {-2^SpAr^(6)AA(r, xr)}] ■ 

(17) 

IV. RESULTS AND DISCUSSION 

We first test our dipole cross section parametrizations 
and vector meson wave functions by comparing them to 



HERA results f5F on diffractive J/^ production that is 
known to be well described by dipole model fits [23J [31] . 
The comparison is quite satisfactory, as can be seen 
from Fig. [2] In addition to the factorized approximation 
(Eq. ([5]), "factorized IPsat" in the figure) that we are us- 
ing in the rest of this paper, also shown is the result with 
the original IPsat parametrization (Eq. ([s]), denoted "IP- 
sat" in the figure). The factorized approximation differs 
from the original one slightly at small Q^, but the differ- 
ence is not significant for our purpose of establishing a 
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FIG. 6: The "nuclear transparency" ratio of cross sections vs. 
xv using the IPsat and IIM parametrizations for = and 
= 10 GeV^ 



FIG. 7: The incoherent cross section integrated over the 
interval 0.1 GeV^ < -t < 0.3 GeV^ divided by the coherent 
cross section integrated over < — t < 0.1 GeV^ as a function 
of Q2 + 



reasonable baseline for computing nuclear effects. 

We note here that the diffractive slope parameters in 
the parametrizations are different, = 4.0 GeV~^ for 
fPsat and = 5.59 GeV~^ for IIM; since these are cor- 
related with the other parameters in the fits leading to 
the parameter values used we do not wish to alter them 
here. Our approximation of a factorized ^-dependence 
with a constant B does not allow us to describe the ob- 
served weak energy and dependence of the diffractive 
slope. The larger B that we use for IIM comes from 
the (To normalization in a fit to inclusive F2 data, and 
also agrees with the observed slopes in inclusive diffrac- 
tion at large /3 and small xp [33] and exclusive p and 
(j) data [34]. The HERA J/*-data, on the other hand, 
has a smaller slope ~ 4 GeV~^ [55]. The t-slope in the 
IPsat parametrization is mostly determined by this J/^*- 
measurement, and an agreement with the larger mea- 
sured slopes for p and (p is obtained by taking into ac- 
count the larger size of the wavefunctions of these lighter 
mesons. 

The differential cross section da'' -^-s- ^^j. ^ _ 
197 (gold) as a function of t is presented in Fig. [3j We 
show the cross section at xp — 0.001 for photoproduc- 
tion. As we performed the nuclear wavefunction aver- 
age leading to Eq. (13) in the approximation where \t\ 
is large, neglecting the coherent contribution, we cannot 
extend our incoherent curves to small \t\. For compari- 
son we show the corresponding "IPnonsat" result where 
the IPsat model is linearized in r'^F(x,r). This curve 
corresponds to the calculation done in Ref. PJ, includ- 
ing both the coherent and incoherent contributions, but 
without the effect of multiple scattering off different nu- 
cleons (i.e. the incoherent cross section is explicitly A 
times the one for a proton). As one can see, the nuclear 
modification due to multiple scattering (resulting mostly 



from the factor e-^-B,AT^(b)[N(r)+N(r')] 
very large. In the full black disk limit of N{r) = 1 this 

and completely supresses the 



factor becomes 



contribution from the center of a large nucleus, leaving 
only an area of « iirdRA ~ A^^^ contributing to the in- 
tegral over by. Thus the cross section in the black disc 
limit behaves as ~ A'^/'^ compared to ~ A in the dilute 
limit, so a large suppression is to be expected. 

We also show in Fig. |3] the coherent cross sections (us- 
ing Eq. (17)). They are also suppressed compared to the 



linearized version (IPnonsat), but not by as much as the 
incoherent one. In the linearized version (as can be seen 
explicitly in Ref. [3] where this case was considered) the 
ratio between the coherent cross section dXt — and the 
incoherent one extrapolated to t = is A. In the IPsat 
model we get 270 (250) and in the IIM model 300 (270) 
at = (Q2 = 10 GeV^). This would make it slightly 
easier to measure the first diffractive dip in the coherent 
cross section, since the background from the incoherent 
process is smaller by a factor of 2 than the linearized 
estimate [1]. 

To demonstrate the nuclear dependence further we 
show in Fig. [4] the ratio of the cross section in a gold nu- 
cleus to that in a nucleoli as a function of . Historically 
this ratio is known as the "nuclear transparency" . Its 
smallness at low energy, similarly to coresponding quan- 
tities in hadron-nucleus scattering, is due to the interac- 
tions of the J l'^ as it propagates through the nucleus. 
The growth of the transparency towards 1 for increasing 
[S] [7] is a demonstration of color transparency (see 
e.g. Ref. [8HI2I [35]). namely that at large the in- 
teracting components of the photon wavefunction are of 
smaller size r and interact weakly. In our framework color 
transparency is automatically present in the fact that 
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FIG. 8: The incoherent cross section integrated over the 
interval 0.1 GeV^ < -t < 0.3 GeV^ divided by the coherent 
cross section integrated over < — f < 0.1 GeV^ as a function 
of xr- 



the dipole cross section approaches zero for r -> 0. In 
Fig. [i] we also show the result (labeled "IPsat, nonsatp") 
of using a nonsaturated dipole-nucleon cross section in 
Eq. ( |13[ ). This corresponds to including unitarity effects 
at the nucleus level but not for a single nucleon. The 
observed nuclear suppression in this unphysical scenario 
is significantly larger than for the saturated full IPsat 
parametrization, showing the sensitivity of the nuclear 
transparency to saturation effects already at the proton 
level. 

The IIM parametrization has a much larger nuclear 
suppression in incoherent diffraction, with the nuclear 
transparency ratio close that of an unsaturated dipole- 
proton cross section. To put this in perspective recall 
that both parametrizations gave an equally good de- 
scription of the elastic cross section measured at HERA 
(Fig.[2|. Since IIM does this with a larger Bp than IPsat, 
we can infer that the typical is smaller, so that the elas- 
tic cross section ~ BpN^ is of the same order. The 
nuclear transparency ratio, on the other hand, depends 
on the total dipole-nucleon cross section ^ BpM ^ a°^/JV 
which is thus larger for IIM. Thus we have a situation 
where both parametrizations have been fitted to inclu- 
sive F2 data^, reproduce well the HERA J/'i' cross sec- 
tion, but differ in their result for incoherent diffraction 
in nuclei. This stresses the importance of performing a 
global analysis of both inclusive and difFractive data to 
constrain the dipole cross sections, and demonstrates the 
utility of eventual incoherent diffractive measurements in 
such an analysis. 



^ Although we have here approximated the original IPsat 
parametrization by factorizing the b-dependence. 
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FIG. 9: The ratio of the coherent (at t = 0) and incoherent 
(at t = —0.5 GeV^, but in our approximation this does not 
depend on t) cross sections to the corresponding ones for a 
proton; normahzed with and A respectively. Plotted as a 
function of A, for xr = 0.001 and = 10 GeV^ 



Figure [5] shows the same Q^-dependence using the 
"gaus-LC" wavefunction. It puts more weight on large 
dipole sizes, leading to a stronger nuclear suppression. 
The cross section ratio typically decreases by ~ 0.04 
from the "boosted Gaussian" wavefunction, but the rela- 
tive structure between the different dipole cross sections 
stays the same. The difference between the cross sections 
themselves is larger, but much of the it cancels in the ra- 
tio. The existing HERA data is not precise enough to 
fully discriminate between different models for the vec- 
tor meson wavefunction, a situation which should also 
improve with planned new DIS experiments. 

The energy dependence of the nuclear suppression 
(again for A = 197) is shown in Fig. [6] for both IPsat 
and IIM parametrizations at = and = 10 GeV^. 
Again we see the larger nuclear suppression in the IIM 
model than in IPsat. The differences in the energy (i.e. 
xp) dependence of the two dipole cross sections are more 
clearly visible in the photoproduction result. This is nat- 
ural, since in the IPsat model the energy dependence 
at the initial scale of the DGLAP evolution (probed at 
smaller Q^) is almost fiat, in stark contrast to the typi- 
cal behavior resulting from BK evolution. At higher 
the difference in the a;-dependence is smaller, although 
there the IPsat-model, driven by the DGLAP evolution, 
turns over to a faster energy dependence. We have not 
extrapolated our curves to higher energies, since there 
is no prospect of experimental measurements. One does 
however see from Fig. [6] that the curves continue to go 
down when extrapolated to smaller ccp. This is to be 
expected since, as discussed previously, one has not yet 
reached the black disk limit. 

In a realistic experimental setup it might be possible 
to detect or veto the nuclear breakup even when the mo- 
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mentum transfer t is not measured very accurately. In 
this case it will be interesting to understand how the rel- 
ative magnitudes of the incoherent and coherent cross 
sections behave as a function of and xp. Generally 
when approaching the black disk limit the coherent cross 
section increases and the incoherent one decreases. The 
relative change shows, however, a smaller dependence on 
and xp than the nucleus/nucleon cross section ra- 
tio. This is shown in our parametrization in Figs. [7] 
andjSj where we plot the the incoherent cross section in- 
tegrated over the interval 0.1 GeV^ < -t < 0.3 GeV^ 
divided by the coherent cross section integrated over 
< -i < 0.1 GeV^ as a function of + M^^^ and 
Xp. Figure [9] further demonstrates the relative similarity 
of the nuclear suppression in the coherent and incoherent 
cross sections. Shown is the A dependence of the ratios 
{da^^^^^/ dt)/{AdaP/ dt) (which, in our approximation, 
is independent of t) and ( dcr^^j^/ dt)/{A'^ da^ / dt)\^^^ for 
Q2 = 10 GeV^ and xp = 0.001. Note that the coherent 
and the incoherent cross sections are normalized by dif- 
ferent powers of A and that width of the coherent peak 
at small t also depends on A. 

Figures [3] and [4] are our main result. Our calcula- 
tion uses as input only well tested parametrizations that 
have been fit to existing HERA data and nuclear geom- 
etry. We work strictly in the small a;-limit which makes 
our formalism simple and transparent. This paper pro- 



vides realistic estimates for the absolute cross sections 
that could be measured in future nuclear DIS experi- 
ments. We have, however, made several simplifying as- 
sumptions in our calculation, the most important being 

a) the factorized impact parameter dependence Eq. ([2|, 

b) the assumption of independent scattering off different 
nucleons Eq. ([T]) and c) neglecting nucleon-nucleon cor- 
relations. Including these effects in a physically correct 
manner and discussing how they could be studied exper- 
imentally is left for future work. As can be seen from the 
values of the nuclear suppression in Figs. |4] and |6j the ef- 
fects of high densities, gluon saturation and unitarity on 
the incoherent cross section are large. Thus incoherent 
diffraction in future nuclear DIS experiments will be a 
sensitive probe of small-x physics. 
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